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We show that in a decay of the form Bd or B s — > P1P27 (where Pi and P2 are pseudoscalar 
mesons), through a flavor changing dipole transition, time dependent oscillations can reveal the 
presence of physics beyond the Standard Model. If Pi and Pi are CP eigenstates {e.g. as in 
Bd — ► -Kstt°7), then to leading order in the effective Hamiltonian, the oscillation is independent 
of the resonance structure. Thus data from resonances as well as from nonresonant decays can be 
included. This may significantly enhance the sensitivity to new physics of the method. If Pi is a 
charged particle, and P2 its anti-particle {e.g. as in Bd — ► ■K + n~ r y), one has the additional advantage 
that both the magnitude and the weak phase of any new physics contribution can be determined 
from a study of the angular distribution. These signals offer excellent ways to detect new physics 
because they are suppressed in the Standard Model. We also show that the potential contamination 
of these signals originating from the Standard Model annihilation diagram gives rise to photons with, 
to a very good approximation, the same helicity as the dominant penguin graph and thus causes no 
serious difficulty. The formalism which applies to the case where Pi and P2 are C eigenstates also 
further generalizes to the case of final states containing multiple C eigenstates and a photon. This 
suggests several additional channels to search for new physics, such as Ksrj'{ri)j, (f>Ks^f etc. We 
also emphasize that the contribution of non-dipole interactions can be monitored by the dependence 
of the mixing-induced CP asymmetry of non-resonant modes on the Dalitz variables. Furthermore, 
using a number of different final states can also provide important information on the contribution 
from non-dipole effects. 

PACS numbers: 11.30.Er, 12.60.Cn, 13.25.Hw, 13.40Hq 

I. INTRODUCTION 



The oscillation of neutral ground state mesons has proved to be a sensitive probe of CP violation and thus a sensitive 
probe for physics at an energy scale well beyond the mass of the meson itself. In recent years the oscillation of the 
neutral Bd meson produced at B factories or hadronic B experiments has provided a means to test the hypothesis that 
the Standard Model (SM) generates CP violation [1] through the Cabibbo-Kobayashi-Maskawa (CKM) mechanism [2]. 
In this approach new physics (NP) becomes evident if the CKM interpretation cannot consistently explain the results. 

Within the realm of B physics, radiative decays resulting from the quark level transitions b — ► 57, where q = d or 
s, have long been recognized [3] as very good channels to look for NP. Indeed the experimental effort [4] to measure 
the rate of b — > 57 both in exclusive and inclusive channels has reached the point where the comparison with the SM 
is dominated by theoretical errors. Further reduction in the theory errors appears rather difficult. Clearly then it is 
advantageous to use in addition to the rates other observables that can reveal new physics. 

One such well known observable is the direct CP partial rate asymmetry. The predicted SM asymmetry [5] for 
Bd — > X s ^f is rather small, about 0.5% and the experimental bound [4] is about a factor of ten above that. For 
Bd — > Xdj the predicted asymmetry in the SM is actually quite large (around 15%) [5] but its branching ratio is 
rather small and also it is experimentally more challenging due to higher backgrounds. Nonetheless, this is clearly a 
promising technique in which efforts are currently being directed. 

Another approach was suggested in [6] (AGS), and has been recently implemented experimentally [7], where the 
oscillations of the neutral B mesons are exploited. Such oscillations in the radiative decays of neutral B mesons 
are suppressed in the SM so if a significant signal of this sort is detected then NP is directly established. Indeed 
the cleanliness of the NP signal in the AGS technique provides one of the important motivations for a "super B 
factory" [8]. 
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The key is the fact that the photon is polarized; the short distance Hamiltonian of the SM makes the very special 
prediction that photons from b (b) decay are dominantly left (right) polarized with the same weak phase. Thus any 
method such as AGS which measures an observable that vanishes with the SM polarized photons is a null test for new 
physics with little dependence on theoretical uncertainties concerning the hadronization of the final state. Recently 
it was pointed out [9] that decays of the form B° — ► P°P°7, where P° is any spin-0 CP eigenstate, can be used to 
probe AGS oscillations in b —* dj transitions. 

The problem with photons from B decay is, of course, that their polarization is not easily measured with current 
detector technology and so indirect methods must be used to probe polarization dependent observables as is the case 
in AGS. Besides AGS some other methods have been suggested to learn about the photon polarization. Consider, for 
example in the case of b — * 57, mesonic decays of the form Bd — > X s -f. Because the initial state has total angular 
momentum J = the hclicity of the photon must be the same as the net hclicity of the X s . In [10, 11] (GP) the 
photon polarization is probed by considering the interferences of various Knir resonances. This approach has an 
advantage that with the four body final state, parity odd observables may be constructed but has a disadvantage that 
the interpretation of the angular distributions requires some understanding of the 1", 1 + and 2 + kaonic resonances. 
Another mode which may be useful to measure the photon polarization in b — > s-f transitions is the recently discovered 
B —* <j>Kj [12]. The photon polarization may also be studied in At decays [11, 13]. 

Another approach suggested in [14, 15] is to "resolve" the photon to an e + e~ pair. This may either be accomplished 
via interaction of the photon with matter through a Bethe-Heitler conversion or internally where the photon is virtual. 
Measuring photon polarization through the Bethe-Heitler conversion may prove experimentally challenging but would 
be a method of general utility. Furthermore, the detection of such e + e~ pairs in the vertex detector at a B factory 
experiment would enable determination of B meson decay positions in decay modes containing only photons in the 
final state. However, since the amount of material in the inner detector must be kept to a minimum in order to 
improve the experimental resolution, the rate of such conversion is rather low. Internal conversion would be present 
and dominate over short distance e + e~ production at low electron-positron invariant mass. The event rate would, of 
course, be smaller than that for direct photons. An additional feature of this channel would be the ability to study 
the CP properties of the short distance electron positron pairs produced at a larger invariant mass [16]. 

In this paper, we will consider a variation on the processes considered in [6] where we will assume that the final 
state consists of two pseudoscalar mesons and a hard photon. For oscillations to occur, the mesons are required to 
be eigenstates of charge conjugation; since they have definite parity, they are also CP eigenstates. The use of a hard 
photon offers the dual advantage that (1) CP violation due to SM bremsstrahlung (discussed in [17]) is suppressed 
and (2) experimental backgrounds are suppressed at higher photon energy. 

Of course, a final state such as Ksn°j is less general than the Kirn 7 states considered in [10], however, in this case a 
simple analysis can lead to powerful conclusions regarding new physics. Indeed, in the case of Bd — > Kstt°j, an initial 
study may be carried out as an extension of the existing analyses for the decay Bd — > K*°j, K*° — > Kstt [7], as our 
discussion shows that all Bd — ► Ksn°-f decays, not only those produced via the K*° resonance, may be included. By 
using such a generalized final state, not only it is possible to increase the statistics, but one can also extract useful 
information regarding the possible contributions to H e g from operators other than the Standard Model dipole. A 
very important characteristic of the mixing-induced CP asymmetry of these modes is that the contribution which 
originates from the dipole term is independent of the Dalitz variables. Thus the key advantage of this class of final 
states over other probes of polarization observables in b — > qj is that in this case the interpretation is relatively clean. 
In the case of b — ► 37, the SM contribution to CP violating observables is only a few percent, so a large signal would 
be an unmistakable sign of new physics. For b — > d^y the SM predicts a much smaller time dependent CP asymmetry 
so that this case may be viewed as a powerful null test. 

In section II, we recall some of the basic issues in radiative decays. Corrections to the dominant dipole H e g and 
their signals are also discussed. Section III briefly reviews the work of [6], whose generalization is the main focus of 
this paper. Sections IV and V contain the main body of our discussions on three-body modes. Section VI briefly 
mentions some generalizations and also presents experimental considerations. Section VII discusses the helicity of 
photons from the annihilation contribution; section VIII briefly discusses effect of non-dipole operators and section 
IX contains a brief summary. The possible complication in the analysis due to the presence of a perturbative phase 
is very briefly outlined in a short Appendix. 

II. BASICS OF RADIATIVE DECAYS 

Let us consider radiative decays of the form B — > P7 where F is either a single meson (e.g. K*) or a multi-particle 
state (e.g. rnrK). The decay is governed by two amplitudes: the decay to right and left polarized photons; the same 
is true for the corresponding decay of the B. We can denote these helicity amplitudes as follows: 
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(1) 



Here we adopt a phase convention for ^l,r such that their phases are equal if the parity of the final state is opposite 
to the internal parity of F. 

At short distances, the photons arise from a radiative transition of the b quark. In most models for this process 
arising from the electroweak scale or higher, it is to be expected that the dominant contribution to the &-scale effective 
Hamiltonian is via the dimension five dipole transition operator, which we will initially assume dominates the process: 



eft 



F q L qo» v - 



75, 



+ h.c. 



(2) 



Here F[ (F R ) corresponds to the amplitude for the emission of left (right) handed photons in the b R — > 9l7l 
ipL — > Qr^Yr) decay, i.e. in the B — > Fjl {B — > F^r) decay. We can relate F[ and F R by defining a parameter ip q , 
which is O i^A in the SM: 

\m b ) 



F[ = F q e i4> t cosV> 9 , 



sinV' 9 , 



(3) 



where 4> q L and <j> q R are CP violating phases. Since the strong interaction respects parity and charge conjugation, this 
model implies that the amplitudes are given in terms of a single complex valued form factor /, so that: 



Ml = -i^/(P$F)[Pinter„al], 

Mr = F R f{$ F ), 

Mr - Ff /(C$ F )[C intCT nal], 

Ml = -F q R *f(CP<P F )[(CP) intcmal ]. (4) 

Here <5>f represents the phase space of the final state F. The sign in front of F^ and F R * arises from the negative 
parity of the initial B° or B° state. Also, Paternal, Gntornai, (C-P)internai are the internal P, C, CP eigenvalues of 
all final state particles which are P, C, CP eigenstates. For example, in case of Bd — > Kstt°^l and Bd — ► Ksn°jR 
decays, we obtain 

^internal = Pjf s P w o = (-1)(-1) = +1, 

Qnternal = C Ks C n C 1 = (- 1) (+1) (- 1) = +1, (5) 

where we use J PC = CT~ for K§ ignoring the tiny CP violation effect in the neutral kaon system. The strong phases 
in / arise from the rescattering between mesons in the final state F. 

In the SM, the contributions are predominantly given by penguin diagrams such as that shown in Fig. 1. By CKM 
unitarity the short distance contribution to this has a CKM phase given by the phase of VtbV t * s . This short distance 
contribution yields predominantly left-handed photons with the right handed component suppressed by m s /mb. This 
right-handed component will also have the same weak phase as the left-handed component. 

The long distance (LD) contributions arise from the c and u penguins. They can have a nontrivial rescattering 
phase distinct from those present in the form factor / mentioned above. An example of a quark diagram contributing 
is shown in Fig. 2, which can also be understood as a rescattering of mesons through processes such as that shown in 
Fig. 3. Fig. 4 represents another LD contribution to radiative decays arising from "annihilation" diagrams; we shall 
discuss this particular contribution later in the paper. 

The precise effect of such LD processes is difficult to calculate reliably. In [18] a detailed estimate in the case of 
B — > V-f is given and it is found that the photons from these contributions are still predominantly left-handed. On 
the other hand, the authors of [15] entertain the possibility that the contribution is of the form Ml = —Mr. 
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FIG. 1: A typical SM radiative penguin graph. 
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FIG. 2: An example of a QCD loop correction which generates the absorptive part necessary for direct CP asymmetry, see [5]. 
The cut is indicated by the dashed line. 



In the case of b — > sj such a contribution will be suppressed with respect to the short distance [21]. However, for 
b — > c?7 the magnitude of such a contribution may be quite appreciable. Phcnomcnologicaly, the primary manifestation 
of such "long distance" contributions would be partial rate asymmetries in B — > K*-f or B — > pj. 

In such cases when there is a contribution which has a strong phase, it can result in a contribution to Eqn. 4 which 
although cannot readily be reliably calculated, may be parameterized as follows: 




FIG. 3: A possible long distance rescattering effect due to on-shell D a D contribution, see reference [18]. 
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FIG. 4: The dominant contribution from the annihilation graph, see reference [19]. 
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-^.f(P$ F )(l + 5 L (P<D F ))[F intcrnal ], 

Ff /(C$ F )(1 + 5 R (C* F )) [^internal], 

-F«7(CP$ F )(1 + <5 L (CP$ F ))[(CP) intcrnal ], 



(6) 



where the 5 terms are arbitrary complex form factors. If such contributions are only from the SM and the photon is 
predominantly the same as the short distance SM effects then: 



where fi is the weak phase between the charm and top penguins and A is a complex valued function of <& F . For 
b — > S7 the small direct partial rate asymmetry from the SM may also be eventually detectable [5, 17]. In the case of 
b — > dj the direct CP asymmetry can be quite sizable. 

In general, new physics should assert itself as an additional contribution to F q and <f> in Eqn. 3 which is different 
from the SM and we will assume therefore that new physics is only manifest at short distances. If this were not the 
case, then probably other signals would be more suited for its detection. Since our methods are generalizations of 
those proposed in AGS [6], let us now briefly review that method. 



Before proceeding it is useful to consider the conditions under which the AGS oscillations occur. These follow from 
the general conditions for oscillating signals in neutral mesons applied specifically to the case of radiative decays. 
Thus for there to be oscillations in the decay B — ► the following conditions are necessary: 

(1) Both B and B must decay to the same exclusive final state V7 (e.g. V = K* , p or higher resonances). 

(2) The photons produced in B — * Vy or B — > V-y must be a mixture of right and left handed hclicities. 

Assuming these conditions are met, we recall that in general we can describe the time dependent wave function of 
the B meson (either Bd or B s ) as [22]: 



S L 

S_R 

Sl 



Sr = 



(7) 



III. AGS OSCILLATION 



\B(t)) = g + \B) + ±g-\B), 
\B(t)) = 9 + \B) + ^9-\B), 



(8) 



where 
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9± = ^" <Ml V 



1 ± e -^mt e iArt 



(9) 



In the case of B d decay, we assume that q/p = e^ M , with </>m = — 20i [24] in the SM, and that AI^ is small. In 
this limit we can write the time dependent decay rates of B d and B d to given final state V-y as: 



^B d ^v~((t) = Y Bd ^v lL {t) + T Bd ^vj R (t) oc e [X Vl + Y Vl cos( Am d t) + Z Vl sin(Am d i)] 



^b^v^) = T B^Vy r (*) + r s rf -y 7R (*) « e ~ Tt \ X Vi ~ Yv 1 cos(Am d t) - Z Vl sin(Am d t)] , (10) 



where 



X v - ( = (|X L | 2 + |Mi?| 2 ) + (|AlL| 2 + |A^fl| 2 ), 
*V 7 = (\M L \ 2 + \Mb\ 2 )-(\Ml\ 2 + \M R \ 2 ), 

Z Vl = -2Im(e l(t > M (Ml~M L + M* r Mr)) . (11) 

Here we sum decay rates for the left-handed and right-handed photon helicity states as we do not distinguish between 
the two. We also define CP violation parameters <Sy 7 and Av 1 as [25] 

l* d ^ Vy us T l^^l = Sv 7 sm(Am d t) + A v - t cos(Am d t). (12) 

The parameter 5y 7 represents mixing-induced CP violation, while Avj represents direct CP violation. From Eqn. 10, 
we obtain <Sy 7 = — Zyy/Xv~ ( and Av-y = —Yy 1 /Xy 1 . 

In AGS it was assumed that Vy was a vector-photon (K*j or fry) final state and that the photon emission was 
described by the short distance Hamiltonian of Eqn. 2. Note that to respect the conditions for oscillations, the K* 
state must be in a C eigenstate which is only true if K* —* Kstt or K* — > K^ir (the latter is unlikely to be 
experimentally useful). Assuming that the polarization information is not available, we must content ourselves with 
summing over polarization and, neglecting long distance contributions, thus obtain the time dependent forms: 



^B d ^v 7 (t) oc e~ rt [1 + xvsin2V> 9 sin<^sin(Am d i)] , 

r s^v 7 W a e- rt [l-xvsin2^sin<^sin(Am d i)], (13) 

where 4> q — 4>m + 4>l + 4>r [26], xv denotes the C eigenvalue of V [27] and the superscript q indicates the quark 
produced in the b — > q-f decay. 

This oscillation therefore allows the extraction of the quantity <Sy 7 = — \v sin2i/) 9 sin^ 9 . In the case of B d — > K*j 
where the short distance contribution to the photon is predominantly right handed, ip s is small and so iSjf* 7 is 
consequently small. The observation of a significant Sk'j would therefore indicate the presence of NP. To the extent 
that the long distance contribution to the photon is right handed as suggested by the calculation of [18], then the 
same is also true of B d — > py. 

Let us now consider the generalization to final states with two pseudoscalars, in particular Ksir°y and ir + ir~y. 
Note that the case of K*y with K* — > Kgir is a special case of Kg^ "/ where the two mesons are on the K* resonance 
while p°j is a special case of 7r + 7r~7 with the two mesons on the p resonance. For clarity, let us remark that any 
kaonic resonance of angular momentum J, that produces Kstt , will have P = (— l)' 7 with C = — 1; thus for any such 
resonance \v — — 1 and so the oscillations for all of them will be identical as they are all governed by Eqn. 13. 

In contrast, for a 7r + 7r~ resonance of spin J, P = {— 1) J and (since CP = +) C = (— l)' 7 . It follows then that 
Xv = {— 1) J > so all odd-J and even-J resonances will have opposite signs in Eqn. 13. 

In the next section we will generalize these results to be independent of the resonance the mesons might go through 
and consider the angular distribution of the pseudoscalars. 



IV. B d — > K s -k°i AND B d — > 

In this section we will contrast the nature of AGS oscillations in the cases of B d — > Kstt ^ [28] and B d — > ir + ir~j. 
The discussion here is easily generalized to the case of other (pseudo-)scalar pairs recoiling against the photon. We 
will consider the more general case including direct CP violation in the next section. 
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The key point to realize is that there is a contrast between the symmetry properties of phase space in the case 
of Kstt with that from tt + tt~. To sec how this arises let us designate the K$ to be particle "1" and the ir° to be 
particle "2" in the first case while we designate n + as particle "1" and tt~ as particle "2" in the second case. For 
Ksn°-f then: 



C*(7k,1,2) - *(7*,1,2), 

P$(7?,l,2) = *(7£, 1,2), (14) 



while for 7T+7T 7, 



C$( 7f ,l,2) = $( 7 «,2,1), 

P$(7«,l,2) = $(7-, 1,2). (15) 

That is to say, under the C and P transformations, each point in the Kstt°j phase space is translated into the 
same point (notwithstanding the distinct photon helicity under P and CP), whereas for ir + ir~j the positions of the 
pions are interchanged under C and CP. 

It is now a simple matter to generalize the AGS formalism to the case where the final amplitude is also a function 
of phase space. Consider first the case of Bd — ► Kstt°j, which at the quark level corresponds to b — > S7. The decay 
amplitude consists of two components which do not interfere with each other corresponding to photons with left and 
right handed hclicities. These amplitudes should depend on the Dalitz plot variables which we will denote: 



\2 



7TIJ J 



51 = (PKs +P- 

52 = (PK S +Plf, 

53 = (P^y+P-Ko) 2 , 

z = . (16) 

S3 + S 2 

In particular the amplitude can be expressed as a function of si and z, where Si is the invariant mass squared of the 
KsTr a system, and z is the cosine of the angle between the Bd and n° in the Kgir frame. 

If we assume that the decay Bd — > ifs7r°7 is governed by Eqn. 2 then the amplitude as a function of s\ and z can 
be written from Eqn. 4 as [29]: 

M R (s u z) = F[*f R ( Sl ,z), 

M L (ai,z) = Fgf L (s u z), (17) 

where the subscript on M indicates the helicity of the photon emitted and the superscript indicates the quark produced 
by the reaction (i.e. d or s). For the charge conjugate decay Bd — > i^s 7r °7 we can likewise write: 

Mr( Si ,z) - Fy R ( Sl ,z), 

M L ( Sl ,z) = F[f L ( Sl ,z). (18) 
Since QCD respects both C and P, we expect from Eqn. 4 

f R (si,z) =J R (si,z) = -f L (si,z) = -J L (S!,Z). (19) 

Thus we obtain 

M L ( Sl ,z) = -F[f( Sl ,z), 
M R (s u z) = F R f( Sl ,z), 
M R (s u z) - F[*f( Sl ,z), 

M L ( Sl ,z) = -F R *f( Sl ,z), (20) 
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where we define a universal form factor f(s\,z) = fn(si,z). Note that in this discussion of the amplitudes at a 
fixed point (s\,z) in phase space, the relative angular momentum between K$ and w° does not enter and therefore 
the quantum numbers (e.g. spin) of intermediate kaonic resonances contributing to the process do not effect the 
conclusions [30]. This is somewhat different from the behavior, discussed below, for the case where the two final state 
mesons are antiparticles e.g. Bd — > 7r + 7r~7. 

The time dependent rates for physical Bd and Bd decays, at a point in phase space defined by (s\, z), and summed 
over photon helicity are therefore given by: 

^B d -^K s n°-y(t, s i, z ) e ~ r * [XK s n°j(si,z) + Y Ksv0j (si, z) cos(Am d t) + Z A - sw o 7 (si, z) sin(Am d t)] , 
T B d ^Ks-K a i^^ Sl > z ) a e " [ X k s w°j(si,z) - Y Ksn o 7 (s!,z) cos(Am d t) - Z Ks7r o 7 (si, z) sin(Am d t)] , (21) 
where, neglecting long distance effects, 



X Ks ^( Sl ,z) = 2(F s ) 2 \f( Sl ,z)\ 2 , 
Y Ks ^, 1 (s 1 ,z) = 0, 

Z Ksv o 7 ( Sl ,z) = -2(F s ) 2 |/( Sl ,z)| 2 sin2^sin<f , (22) 



and 4> s = 4>m + 4> s l + 4> s r is the weak phase. 

Thus, for each value of s± and z, the CP asymmetry is: 



r B d ^K s ir«~< ~ ^B d ^K S TT°~t 
T B d ^K s ir<'~< + ^B d ^K S Tt°~t 



= + sin 2ip s sin 4> s sin( Am d t) . (23) 



Note that this expression is true whether the Ksir is produced by the decay of a resonance or is nonresonant. 
Furthermore, the fact that the effective Hamiltonian of Eqn. 2 produces the photons implies that this asymmetry 
does not depend on sj. In effect each point in phase space is a separate oscillation experiment which shows the 
same oscillator behavior given by Eqn. 23. In practice this means that one may add together all events of the 
form Bd — > Kstt°j regardless of whether they are produced at the T^T* (892) resonance, some other resonance (e.g. 
if2(1430)) or from a nonresonant source and determine the single quantity SB d —>K s ir°-y = + sin 2-0 s sin </> s as a result. 
Within the SM, we obtain S Bd _ Ks ^o 1 rts — (2m s /mb) sin 2<j>\ as ip w m s /rrib and s = — 20i + 0(A 2 ), are expected 
where A w 0.22 is the sine of the Cabibbo angle. In terms of the individual parameters this gives the lower bounds: 



|sin2V> s | > |5 Bd ^if sw o 7 | , 
|sin0 s | > \S Bd ^ Ksn o 7 \. (24) 

Deviations from this picture of uniform oscillation as a function of phase space would indicate contributions to the 
radiative decay other than the pure dipole transition of the effective Hamiltonian of Eqn. 2. 

Consider now the case of Bd — ► ir + ir~^. Again for the Bd decays we can define the Dalitz variables: 



and thus write the amplitudes in the form: 



51 = (Pir+ +Ptt-) 2 , 

52 = (Prr+ +Pj) 2 , 

53 = (f> 7 +Pn-) 2 , 

S3 - S 2 

z = ! ' 

S3 + s 2 



(25) 



M L (si,z) = F£g L ( Sl ,z), 

M R (s!,z) = F^g^suz), 

M R (s u z) = Ftg R (s u z), 

M L (s u z) = Ffg L ( 8l ,z). (26) 
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As before g^ = —gn but in this case in applying charge conjugation to get from g^ to g L we need to interchange the 
coordinates of the n + and n~ . Thus s 2 <-> s 3 so z ^> —z, therefore [31]: 

9l(si,z) = -g R (s!,z) = -g(s!,-z), 

z) = -g R (s!,z) = -g{si, z). (27) 

For particular partial waves of the 7r + 7r~ system with angular momentum L, g(si, z) = {— l) L+1 g(s\, — z). In general 
g will be a mixture of even and odd L, so g has no particular symmetry under z <-» — z. Note also that 5 will in 
general have a nontrivial CP-even phase which depends on s\ and z. The time dependent rates at a point in phase 
space defined by (s±, z), and summed over helicity will be given by: 

(s-l, z) sm(Am d t)] , 

r B <J - w + w - 7 (*' *ii z ) « e_r * [^7r+ w - 7 (si) z ) - Y 7,+n- 7 (si,z) cos(Am d t) - Z 7r+7r - 7 (si, z) sin(Am d t)] . (28) 

Here: 

X^-^suz) = (F d ) 2 [\g( Sl ,z)\ 2 + \g( Sl ,-z)\ 2 ], 
Y^-^suz) = (F d f[\g(s u z)\ 2 -\g(s u -z)\ 2 }, 

Z„ + „- 7 ( 8l ,z) = -2(F d ) 2 [Rc(.g( Sl ,-z) 5 *( Sl ,z))sin^ + Im( ff ( Sl ,-z) ff *( Sl ,z))cos^]sin2^, (29) 

and <ft d = (f>M + fit + 4>r is the weak phase; in the SM </>m ~ — 2</>i and = <^ <~ <^>i. 

At each point in phase space, g(si, z) and g(si, —z) are, in general, independent complex numbers. Unlike the case 
of i^s7r°7, this case does, however, allow the possibility of extracting <j> d and ip d separately up to discrete ambiguities. 
This can be achieved without any assumption about the resonant structure of the Tr + n~-f final state, as we now 
demonstrate. Note that our argument requires regions of phase space where tt + tt~ partial waves with different 
angular momentum interfere, otherwise g(s\,z) — ±g(si,—z) everywhere and consequently there is no additional 
information compared to the KsTr "f case. 

For a given value of Si and z let us define: 

« = cos2^, a 

« = sin2^sin^, b = -2 ^^^ ; :$ , (30) 
w = sin2^cos^, c = 



where a 2 + b 2 + c 2 = u 2 + v 2 + w 2 — 1. Recalling that the experimental observables, X n + n - 7 , Y n + n - 1 and Z K+7r - 
are functions of the phase space, we can relate these via: 

V = y = a ' 

C+ = ^ 1 - = bv, 

y y 



c = — 1 = cw. (si) 



Along the line z = 0, we find = 7/ = while £ + = —v. Once z; is known, we can learn w from cases where z ^ 0: 

- - , = , (32) 

2 / J I ^2 

,2 



From here we can determine (/> d and ^ d through 



sin 2 2?/^ = w 2 + w 2 , 

</> d = arg((w; + w)/sin(2V> d )) . (33) 
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There are 8 solutions for ip d . Each of them has two <j) d values corresponding to positive and negative solutions for w. 
Thus there are 16 valid (ip d ,<j> d ) combinations in total. 

In the preceding discussion, we have left the functions gL,R arbitrary. In practice, the tt + it~ amplitude is very likely 
to be dominated by low-lying resonances with well known masses, widths and quantum numbers. This knowledge could 
facilitate a more constrained fit to the data. In case Bd — > p (—* ir + ir~)j is dominant, we expect g(s\, z) = g(si, —z) 
as the 7r + 7r _ system is in the L-odd state. Substituting this into Eqn. 29, we obtain <S Bd ^ p o 7 = + sin 2tp d sin 2<p d . 

Throughout this discussion we have, again, assumed dominance of the dipole operators in H c g, as in Eqn. 2. To 
the extent that this assumption holds, ip and <fr would be independent of the Dalitz variables, s\ and z, as before. 
Conversely, variation of ip and <f> with s\ and z would shed light on the contribution from non-dipole interactions. 

V. B a -> K+K-f AND B a -f Ks^iv^h 

Let us consider now the analogous B s decays. In this case we cannot assume that Ar w 0; then the time dependent 
decay rates to a final state Fj are given by: 

^B s ^Fy(t) oc e~ r * [X Fl cosh(iAr s i) + Y Fj cos(Am s t)+ 

Z Fl sm(Am s t) + W Fl sinh(±AT s t)], 
r s s ^F 7 (*) x e_rt [ X Fi cosh(iAr s t) - Y Fl cos(Am s t)- " 1 > 

Z Fl sm(Am s t) + W Fl sinh(iAr s t)] , 

where 



X Fl = (\Ml\ 2 + \Mr\ 2 ) + (\M l \ 2 + \M r \ 2 ), 

Y Fl = (\Ml\ 2 + \Mr\ 2 )-(\M l \ 2 + \M r \ 2 ), 

Z F - ( = -2Im(e^ M (M* l M l +M r Mr)) , 

W Fl = -2Rc(e^ M (MIM L + M* R M R )) ■ (35) 

Here </>m ~ 0(\ 2 ) is expected for the B s -B s mixing in the SM (in the usual phase convention). 

Let us first consider the decay to two mesons which are self-conjugate such as B s — > Ksir°~f. The expression Eqn. 34 
also describes the time dependent decay rates for physical B s and B s decays at a point in phase space defined by 
(si, z), and summed over photon helicity. Neglecting long distance effects, the factors X, Y, Z and W are given by 



X Ksw o( Sl ,z) = 2(F d ) 2 \h( Sl ,z)\ 2 , 

yK s n<>{si,z) = 0, 

Z Ks7r °(si,z) = -2(F d ) 2 \h(s 1 ,z)\ 2 sin2iP d sm<j) d , 

W KsV o(s u z) = -2(F d ) 2 \h( Sl ,z)\ 2 sm2^ d cos(f> d . (36) 

In this case, if the value of W can be isolated, both angles 4> d and ip d can be determined, up to discrete ambiguities. Of 
course if AT/T is relatively small, this may prove difficult in practice. As in the case of Bd — + Kstt°j, the asymmetries 
are independent of Si and z. 

In the case of B s — > K + K~j, the phase space dependent oscillation is again given by the expressions in Eqn. 34. 
The X, Y, Z and W terms are given by: 



X K+K -(s u z) = (F s ) 2 (\d( Sl ,z)\ 2 + \d( Sl ,-z)\ 2 ), 

Y K+K -(s u z) = (F s ) 2 (\d( Sl ,z)\ 2 -\d( Sl ,-z)\ 2 ), 

Zr+ k- (si,z) = -2(F S ) 2 {Re (d(si,-z)d*(si,z)) sin 4> s + Im (d(si, -z)d*(si, z)) cos <f> s } sin 2ip s , 

W K+K -(si,z) = -2(F s ) 2 {Re(d(si,-z)rf*(si,z))cos ( / ) s -Im(d(si,-z)d*(si,z))sin0 s }sin2V' s . (37) 

From the X, Y and Z terms one can determine 4> s and tp s as described above in the Bd tt + it~"/ case. In addition, 
if the W term can be determined it allows for the resolution of the ambiguity between <j> s , ip s and <fi s , —tp s . 
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its it's 7 


7T + 7T 7 






Bd/B d 


b — > s-y 


6 — ► d7 


6 — > d7 


b — ► d7 


6 — > d7 


B a /B s 


b — > c?7 


6 — > S7 


6 — > S7 


b — > S7 


6 — ► S7 



TABLE I: Final states which can be used to probe b — > S7 and & — > d"f transitions in Bd and B 3 decays. This list is not 
exhaustive; in particular other neutral (pseudo-) scalar particles (r/, rj ', /o) may be used in the place of 7r°. 

VI. ANALOGOUS CASES, GENERALIZATIONS AND EXPERIMENTAL CONSIDERATIONS 

In the above discussion we have considered some instances of Bd — > P1P27 and B 8 — » P1P27 where Pi represent 
scalar or pseudoscalar states. The mode of analysis depends on the nature of P1P2. 

If there is no relation between Pi and P2 and they are not self-conjugate then there will be no oscillations. For 
instance in the case of Bd — ► K + tt~j we can tell from the final state whether the initial state is B or B so no quantum 
interference is possible. Final states of this sort do however give the simplest way to determine if there is direct CP 
violation at the quark level. 

If Pi and P 2 are both eigenstates of charge conjugation, then the mode of analysis is the same as Bd — ► Ksn°j. In 
this case we can learn the product sin2i/) 9 shi(/> 9 , where the superscript q represents the quark produced in b — > q-f. 
For B s oscillations, if the W term can be extracted then we can learn <fi q and ip q separately. In all cases, the data can 
be integrated over si and z. Some other final states of this type, relevant to Bd decays, are Ksff'^, KsfoJ, n°rj^j, 
7r°/o7, For most of these modes there is currently no experimental information; they are therefore in urgent need of 
investigation. Cases that may be of special interest are B° Ksr/'j and B° — ► Ksrjj. Comparisons to the pattern 
of branching fractions in two body hadronic B decay suggest that the nonresonant contribution to B — > Kr/'^y ought 
to be larger than that for either B — ► if 777 or B — > Kir°j. Note also that B — > if 777 has recently been observed with 
branching ratio w 7 x 10 -6 [32]. Therefore it is possible that an appreciable data sample for if 57/7 may already be 
available with the current B factory statistics. 

Note that there is a set of special cases of these modes, where Pi = P 2 , for instance, if s if 57 and 7r°7r°7 (the latter 
unlikely to be of practical use unless a very high luminosity allows us to use n° — ► e + e~-f decays that provide vertex 
information). In these cases, Bose-Einstein statistics further constrain the P1P2 system [9]. 

As explained in Sec. IV, the relative angular momentum between Pi and P2 does not enter; only the intrinsic charge 
conjugation quantum numbers of Pi and P2 affect CP asymmetries. This is also valid in more general cases with 
more than two particles, e.g. B — > P1P2P37, where Pi, P2 and P3 are eigenstates of charge conjugation. Specific 
examples of this type are B — > ifg7T°?T 7, ifgifsifs7, ifsifs?/'^, and so on. 

Indeed, as long as all the final state particles are eigenstates of charge conjugation, regardless of their spin or other 
quantum numbers, Eqn. 23 applies since charge conjugation does not map one point of phase space to another. A 
case of particular interest is B° — > Ks<jry. This final state has the practical advantage that the <f> can be observed 
as K + K~ , allowing a simple determination of the B decay vertex, rather than the extrapolation which is necessary 
when using if 5 — > 7T + 7r~. Note that there is an indication that this mode has a branching fraction of <~ 5 x 10 -6 [12]. 
Other analogous cases of interest are Ksfq and Ksowy. 

If Pi and P2 are not self-conjugate but are anti-particles of each other then the mode of analysis is the same as 
7r + 7r~7 above. In this case by considering the time dependent Dalitz plot we can separately determine <f> q and ?p q - 
Final states of this sort are 7r + 7r~7 and K + K~-f. Table I shows the various final states which may result from Bd 
and B s decay and which quark-level process they are sensitive to. 

Let us now consider a few of these cases which are likely to be of greatest experimental interest. In the case 
of Bd decay preliminary studies have been done of the AGS mode Kstt°j on the K* resonance [7]. Part of the 
importance of the discussion in this paper is that this data may be combined with nonresonant Ksn ^ decays and 
also with signal from other resonances. Since we are interested in the underlying two-body process b — > qj, we assume 
that experimental cuts will be imposed to discriminate against bremsstrahlung and other possible background. At 
a P-factory experiment, these cuts will typically include a requirement on the center-of-mass energy of the photon; 
since the B meson is almost at rest in the T(45) rest frame this is equivalent to requiring a hard photon and will 
remove most bremsstrahlung events. Further cuts may include other requirements on the P1P27 phase space. In 
addition to reducing the dominant experimental background from continuum e + e~ — > qq, q = u,d, s,c events, these 
requirements can remove the background caused by two body hadronic B decays followed by radiative hadronic decays. 
For example, the decay chain Bd — ► cjKs, to — ► 7r°7, contributes to the final state Ks^-f. This background decay has 
a product branching fraction of 2.5 x 10~ 7 [23], which is small, but not entirely negligible compared to that for the 
signal. However, this background can be removed with a requirement on the invariant mass of the 7r°7 combination. 
Similar backgrounds should be considered for each final state. 

Replacing the n° with an jj(') is more challenging experimentally, although the decay B — > if 777 has recently been 
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observed [32]. These decays should measure the same quantity if the decay is controlled by Eqn. 2. It would however 
be of some importance to use these modes as a check that new physics proceeds through this dipole operator. 

In the Bd system K$Ks7 monitors sin cj) d sin 2\p d for the b — > d-f transition. Of course this case may also be subject 
to a significant amount of direct CP violation at the quark level which can also be measured in the usual way in a time 
(in)depcndent analysis, and could be further checked via charged B decays such as B + — ► K + Ksj or B + — > 7r + 7r°7. 

The P + P~-f states Bd — > n + ir~-f and B<j — > K + K~j could in principle be used to separately obtain (j) d and tp d . 
The — ► K + K~j case is, however, OZI suppressed and so is unlikely to have a significant branching ratio. The 
same is true for B s — ► 7r + 7r~7. Experimentally, good A"/7r separation is needed to isolate Bd — > 7r + 7r~7 from the 
decay Bd — > A" + 7r~7 which may have an order of magnitude greater branching ratio. 

Observing £? s oscillations at hadronic £> machines will of course be challenging. The simplest of the modes in 
Table I to measure is probably B s — ► K + K~-f which again is sensitive to the b — > S7 transition. 

VII. ANNIHILATION CONTRIBUTION 

The key feature of oscillations in B q — ► B 7 is that they only may take place if ?/> 7^ 0. In the SM, this reaction 
takes place through the penguin in Fig. 1. The photon due to this penguin process is exactly left handed in the limit 
of massless q; since in that case q would be left handed, so conservation of helicity would imply that the the photon 
would likewise necessarily be left handed. More generally, in the SM this process gives rise to tp w m q /mb. 

Another source of right handed photons which could therefore potentially produce a signal is the annihilation 
diagram shown in Fig. 4. This process would only contribute as a background to be present in b — > dj processes in 
Bd decays and to b — > sj processes in B s decays. These processes are subject to large nonperturbative effects when 
the photon is radiated colinearly with the initial light quark. This enhancement, however, is only in the left handed 
photon channel. To see this consider the amplitude of the relevant annihilation graph (d) where we have applied a 
Fierz transformation to the W propagator: 

M ann ^ 9 o 

— (dWi - 1> d )l»Lb) (d^Lu) , (38) 

Pd~ ^1 ' P d 

where pd is the momentum of the initial d quark, L = (1 — 7 5 )/2 and R = (1 + 7 5 )/2. We can rewrite this in the 
following form: 

Man„ cx 9 E„q v R] 7 "6) (d 7 „Lu) + -3-^ (dWrfLb) (d^Lu) . (39) 

Pd ^1 ' Pd Pd ' Pd 

where fs is the B decay constant and is the polarization 4-vector of the photon. The first term is by itself gauge 
invariant. The factor in square brackets is the dipole operator for the emission of left handed photons. This term 
is enhanced by the propagator in front because this quantity in the rest frame of the B meson is approximately 
l/(E^Ed) where Ed is expected to be small w Xqcd- As discussed in [19] this term corresponds to emission of a 
colincar photon by a light (initial) quark and because of its singular nature in perturbation theory it is expected to 
receive non-perturbative corrections. Since it has the same photon helicity as the penguin produced photons but a 
different CKM phase, it will alter the magnitude and phase of F£ but will not contribute to and therefore will not 
much affect V' 9 - 

Actually, the singular nature of the first term in Eqn. 39 is a consequence of the very simple non-relativistic model 
used for the bound state. Indeed in more sophisticated discussions, that quantity is not singular anymore but has 
a well defined light-cone expansion in 1/E 1 [20]; the important point relevant to this paper is that the left-handed 
nature of the simplistic picture above survives the improved theoretical treatment. 

The second term is not gauge invariant but must form a gauge invariant set when added to the other bremsstrahlung 
diagrams. This term is proportional to the 4-momentum of the light quark and so the energy of the light quark in the 
numerator will tend to cancel the denominator. Therefore, though these graphs will produce photons of both helicity 
they are not enhanced by a light quark propagator and so are expected to make only a negligible contribution. 

Thus we conclude that the annihilation contribution does not cause any serious difficulty to the analysis above. 

VIII. NON-DIPOLE OPERATORS 

A potential complication for the use of this method as a test for new physics is the SM contribution to radiative 
decays through non-dipole operators [33], for instance, the penguin b — > sjg. Such processes generally do not fix the 
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helicity of the photon and so would result in a SM contribution to S. Also the precise contribution of such processes 
to particular final states is difficult to calculate reliably. Fortunately, the experiments being suggested here can give 
a handle on the extent of these non-dipole contributions and the presence of the latter need not represent a serious 
limitation to the application of our method to search of NP. 

The different operator structure in H c g would mean, that in contrast to the pure dipole case, S would depend on 
the kinematics and composition of the final state. For example, for the case of modes such as Ksn ^,!]')^ (see the 
discussion after Eqns. 23 and 24) the presence of non-dipole contributions would, in general, make the asymmetry S 
depend on the Dalitz variables si and z. Thus, the contamination from non-dipole terms in H e g may be estimated by 
fitting the experimental data on S to a suitable parametrization of the dependence on si and/or z. A difference in the 
values of S for two resonances of identical J PC would also indicate non-dipole contributions. Similarly for the n + ir~-f 
mode, ip and 4> would depend on s\ and z. For a better understanding of these non-leading effects, it is also useful to 
measure S for a number of different final states and kinematic ranges and see if there is variation. For instance, in 
the case of b — > some decays which should be studied are those suggested above of the form B° — > jKs + neutrals 
where neutrals could include 7r°, rf'\ p°, lo, 0, Kg, ... or any combination of such mesons. 

IX. SUMMARY 

We have extended the work of [6] so that data from other resonances as well as nonresonant contributions to final 
states such as ifs7r°7, Ksfo^f, K + K~j, tt + tt~j, Kgr]^ '7 and KgKg^ can be included for the b — > s or b — > d 
transitions relevant to Bd or B s decays. This should significantly improve the effectiveness of testing the Standard 
Model with mixing- induced CP violation in radiative B decays. Indeed, this not only helps in reducing the statistical 
errors, since no separation between resonant and nonresonant events is needed, it should also help in reducing the 
systematic errors. For the Kstt ^ type of final state an improved determination of the product sin 2ip q sin 4> q , where 
(f) q is the weak phase and ip q monitors the photon helicity, is possible. For final states such as ir + ir~j and K + K~-f, 
separate determination of each of these two quantities is possible. 

A key feature of these modes is that the dipole interaction in H e g gives rise to a mixing-induced CP asymmetry 
which is independent of the Dalitz variables. Recall that in many models of New Physics it is the (dimension five) 
dipole interaction that is likely to be dominant. Thus the study of the CP asymmetries of these modes should be 
useful in searching for New Physics even in the presence of non-dipole SM contributions. 

In passing, we briefly recall the hierarchy of CP asymmetries in radiative B decays expected in the Standard Model. 
Assuming the dipole Hamiltonian dominates, for b — > s the mixing-induced CP asymmetry is expected to be 0(3%) 
and the direct CP asymmetry [5] should be around 0.6%. For b — > d, the direct CP asymmetry is expected to be 
around 15% whereas the mixing-induced CP asymmetry is « 0.1%, making it into a very interesting (essentially) null 
test of the SM. 
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Appendix: The Effect of a Perturbative Phase 

As discussed in [5] the b — > dj transition is expected to receive contributions with different weak and strong phases 
which could lead to appreciable direct CP violation at the quark level. If we take this effect into account in the case 
of Bd — > 7r + 7r~7, the modified values of A, Y and Z are: 

A(7r+7r- 7 ) = F 2 {\\g( Sl ,z)\ 2 + \g( Sl ,-z)\ 2 } (1 + S 2 ) + 25 cos 2 V [cos( M - v)g{si,z) + cos( M + v)g( Sl , -z)}} , 
y(7r+7r-7) = F 2 { [\g{ Sll z)\ 2 - \g( Sl , -z)\ 2 \ (1 + 28 cos 2 V cos^cos^)- 

[\g(si, z)\ 2 + \g(si, -z)\ 2 ] 25 cos 2 V>sin^sin^} , (40) 
Z(TT + Tr-j) = -2(F d ) 2 {Re(g(s 1 ,-z)g*(s 1 ,z)) [sin^ + 5 cos p sin(<j) d + v)} + 

Im (g(si, — z)g*(s\, z)) [cos</> d <5cos/zcos(</> d + v)\ } s\n2ij) d . 
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